Cosmological equivalence principle and the weak— field limit 
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The strong equivalence principle is extended in application to averaged dynamical fields in cosmol- 
ogy to include the role of the average density in the determination of inertial frames. The resulting 
cosmological equivalence principle is applied to the problem of synchronisation of clocks in the ob- 
served universe. Once density perturbations grow to give density contrasts of order one on scales 
of tens of megaparsecs, the integrated deceleration of the local background regions of voids relative 
to galax;ies must be accounted for in the relative synchronisation of clocks of ideal observers who 
measure an isotropic cosmic microwave background. The relative deceleration of the background 
can be expected to represent a scale in which weak-field Newtonian dynamics should be modified 
to account for dynamical gradients in the Ricci scalar curvature of space. This acceleration scale is 
estimated using the best-fit nonlinear bubble model of the universe with backreaction. At redshifts 
z < 0.25 the scale is found to coincide with the empirical acceleration scale of modified Newtonian 
dynamics. At larger redshifts the scale varies in a manner which is likely to be important for under- 
standing dynamics of galaxy clusters, and structure formation. Although the relative deceleration, 
typically of order 10~^°ms~^, is small, when integrated over the lifetime of the universe it amounts 
to an accumulated relative difference of 38% in the rate of average clocks in galaxies as compared to 
volume-average clocks in the emptiness of voids. A number of foundational aspects of the cosmo- 
logical equivalence principle are also discussed, including its relation to Mach's principle, the Weyl 
curvature hypothesis and the initial conditions of the universe. 
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I. INTRODUCTION 

The strong equivalence principle (SEP) stands at the 
conceptual core of general relativity, as a physical prin- 
ciple which limits the choice of our physical theory of 
gravitation among all possible metric theories of gravita- 
tion one can construct. In this paper I will argue that 
the ramifications of this principle have not been fully 
explored, and that its physical interpretation requires 
further clarification to deal with the dynamical proper- 
tics of spacetime inherent in Einstein's theory when the 
nonequilibrium situation is considered. In particular, the 
problem of how to synchronise clocks in the absence of a 
spacetime background with specific symmetries does not 
have a general solution in general relativity. In this paper 
I will show that at least for universes which began with 
a great deal of symmetry, as ours did, the reasoning of 
the equivalence principle can be extended: the average 
regional density provides a clock in expanding regions. 
This has particular consequences for cosmological models 
and the definition of gravitational energy. It underlies the 
author's proposal that dark energy is a misidentification 
of cosmological gravitational energy gradients in an in- 
homogeneous void-dominated universe [J, [2i 9 ■ Broader 
foundational consequences may also follow. 

To set the scene, itpays to recall that historically the 
equivalence principle^, and indeed general relativityf^. 
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was formulated at a time before the dynamical properties 
of spacetime were understood. The conceptual route that 
Einstein took began with the weak equivalence principle 
or the principle of uniqueness of free fall, known since 
the experiments of Galileo, that all bodies (subject to 
no forces other than gravity) will follow the same paths 
given the same initial positions and velocities. Realising 
that this observational statement embodies a feature of 
universality of the gravitational interaction, Einstein cre- 
ated a theory in which gravity is a property of spacetime 
itself. 

Einstein's identification of what the true gravitation 
field should be began 101 years ago with first identifying 
what it is not, based on thought experiments with ele- 
vators, concerning what motions of particles cannot be 
distinguished operationally. His 1907 principle of equiva- 
lence [J| may be translated as follows: All motions in an 
external static homogeneous gravitational field are iden- 
tical to those in no gravitational field if referred to a uni- 
formly accelerated coordinate system. A uniformly ac- 
celerated reference frame may be found operationally in 
empty Minkowski spacetime by firing rockets; if matter is 
the source of the true gravitational field then such choices 
of frame cannot represent gravity. 

More generally, since special relativity with nongrav- 
itational forces appears to always be valid in small re- 
gions, one should always be able to get rid of gravity 
near a point. This is embodied in the SEP: At any event, 
always and everywhere, it is possible to choose a local in- 
ertial frame such that in a sufficiently small spacetime 
neighbourhood all nongravitational laws of nature take on 
their familiar forms appropriate to the absence of grav- 



2 



ity, namely the laws of special relativity. This means 
that gravity is made to be universal, as it is contained in 
spacetime structure. The true gravitational field strength 
is encoded in the Riemann curvature tensor, and is deter- 
mined regionally by the tidal effects of geodesic deviation. 

One of the most profound and difficult consequences 
of the SEP is that gravitational energy and momentum 
cannot be described by a local density, and so are not 
local quantities. General relativity overcomes the nonlo- 
cality problem of Newtonian gravity: there is no action at 
a distance. General relativity is an entirely local theory 
in the sense of propagation of the gravitational interac- 
tion, which is causal. However, the background on which 
the interaction propagates may contain its own energy 
and momentum, when integrated over sufficiently large 
regions, and this has to be understood in the calibration 
of local rods and clocks at widely separated events. 

Whereas the calibration of rods and clocks is math- 
ematically determined by invariants of the local metric, 
and the spacetime connection which relates local invari- 
ants at widely separated events, in practice we cannot 
analytically solve Einstein's equations for the most gen- 
eral distribution of matter to unambiguously determine 
the metric and its connection. A slicing of spacetime into 
hypersurfaces, and a threading of these hypersurfaces by 
timelike worldlines of observers or by null geodesies, is 
inevitable for any operational description of spacetime in 
terms of rods and clocks. Such splittings of space and 
time, together with additional symmetry assumptions, 
are also necessary for analytically solving Einstein's equa- 
tions in particular cases, or more generally for numerical 
modelling. 

The definition of quasilocal gravitational energy and 
momentum ^ then turns out to depend on the choice of 
slicing, associated surfaces of integration, and the iden- 
tification of observers that thread the slices. These pro- 
cedures are inherently noncovariant and nonuniquc, and 
many questions of naturalness of any particular defini- 
tion inevitably arise. (See Ref. [2] for a recent discus- 
sion.) We have the dilemma that a spacetime split in- 
evitably breaks any given particle motion into a motion 
of the background and a motion with respect to the back- 
ground; and this may involve a degree of arbitrariness. 
The viewpoint that will be adopted here is that since 
quasilocal gravitational energy gradients have their ori- 
gin in the equivalence principle, the primary criterion 
for making canonical identifications of physically relevant 
classes of observer frames is that the equivalence principle 
itself must be properly formulated, and respected, when 
making macroscopic cosmological averages. 

If the SEP is applied to macroscopic objects then 
strictly speaking we can only apply it to systems in which 
the gravitational interaction is ignored. Yet we implic- 
itly apply the SEP to scales at least as large as galaxies 
which are treated as particles of dust in an expanding 
fiuid - with an expansion rate given by a Friedmann- 
Lemaitre-Robertson- Walker (FLRW) model - subject 
to additional Newtonian gravitational interactions. The 



Newtonian approximation, which is made in present-day 
numerical simulations of structure formation must, by 
the rules of general relativity, presuppose that a suffi- 
ciently large static Minkowski frame can be found. 

I shall argue in this paper that to correctly derive a 
Newtonian limit, without prior assumptions about the 
cosmological background geometry, we must first cor- 
rectly apply the SEP. If galaxies are to be treated as 
particles of dust we must address the following ques- 
tion: given that the background is not static, what is 
the largest scale on which the SEP can be applied? An 
attempt to answer this question, which is unavoidable if 
we are to consistently apply the principles of general rel- 
ativity to cosmology, means that we have to deal with 
the relationship of inertial frames to averages of mat- 
ter fields and motions. In other words we must address 
Mach's principle, which may be stated as [1, [§| follows: 
"Local inertial frames are determined through the distri- 
butions of energy and momentum in the universe by some 
weighted average of the apparent motions". The leading 
questions in this statement are what is to be understood 
by "local" , and what is the "suitable weighted average" ? 

The problem with any process of averaging in general 
relativity is that by coarse graining we can lose informa- 
tion about the calibration of local rods and clocks within 
the coarse grained cells relative to average quantities. 
Rods and clocks are related to invariants of the local met- 
ric, can vary greatly within an averaging cell, and will not 
in general coincide with some average calibration of rods 
and clocks once averaging cells become sufficiently large 
in the nonlinear regime of general relativity. Galaxies, for 
example, contain supermassive black holes, in whose local 
neighbourhood the determination of proper lengths and 
times for typical observers differs extremely from typical 
observers in the outskirts of a galaxy. 

It is commonly believed that as long as we are "in the 
weak-field limit" we do not have to worry about compli- 
cations such as the extreme ones posed by black holes. 
However, the weak-field limit is always taken about a 
background, and once inhomogeneities develop in the uni- 
verse there are no exact symmetries to describe the back- 
ground. One set of uniformly calibrated rods and clocks 
is no longer sufficient to describe the background itself. 
Given an initially homogeneous isotropic universe with 
small scale-invariant density perturbations, once one is 
within the nonlinear regime of structure formation, ho- 
mogeneity and isotropy are only defined in a statistically 
average sense. 

Within a cell of statistical homogeneity - which can be 
taken to be [Ij of the same order as the baryon acous- 
tic oscillation (BAO) scale, 100/i~^Mpc, h being the di- 
mensionless parameter related to the Hubble constant 
by — 100ft. km sec^^ Mpc~^ - there are density con- 
trasts of order unity over scales of tens of megaparsecs. 
Since the universe is inhomogeneous over these scales not 
every observer is the same average observer. Different 
classes of canonical observers are required to interpret 
cosmological parameters 1]. In particular, we and the 
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objects we see are typically in galaxies in locally non- 
expanding regions which formed from density perturba- 
tions which were greater than the critical density. In 
such regions local spatial curvature can differ markedly 
from the volume-average location in voids, where space 
is freely expanding. Small differences of spatial curvature 
and rates of clocks can accumulate to give large differ- 
ences over the lifetime of the universe [l!]. Dynamical 
gradients in the curvature of space are a physical reality 
which must be properly understood. 

In this paper I will estimate the scale of the small rela- 
tive deceleration of the background in regions of different 
density, by proposing an extension of the SEP to cosmol- 
ogy. I will demand an equivalence between the particu- 
lar example of geodesic flow of a congruence of particles 
"at rest" in a dynamically expanding universe, whose av- 
erage volume expansion is governed by an average over 
all masses and motions of the particles, and the equiva- 
lent "volume-expanding" motion of point particles in a 
Minkowski space. In other words, the local Ricci scalar 
curvature due to the volume average of pressureless dust 
can always be "renormalised away" on sufficiently small 
scales, but in a way which may lead to relative recalibra- 
tions of local rods and clocks between different spacetime 
regions. 

Equivalently, for volume expansion we cannot locally 
distinguish between particles at rest in a dynamic space- 
time, and particles moving in a static spacetime: the two 
situations are equivalent in a sense which deserves the 
designation cosmological principle of equivalence. This 
might be viewed as a further Machian style refinement 
of the notion of inertia. Although we measure geodesic 
deviation, in terms of the scalar curvature part of the 
Riemann tensor and volume expansion, we are unable to 
distinguish whether the geodesies are deviating because 
of local accelerations of particles in a static space, or 
whether the particles are "at rest" in an expanding space 
which is decelerating due to the gravitational attraction 
of the average density of matter. 

Historically, it might be said that although Einstein 
was conceptually guided by Mach's principle, he never 
quite succeeded in fully implementing it in general rel- 
ativity, because when he first formulated the theory he 
did not fully appreciate the dynamical nature of space- 
time. His first attempt to study cosmology indicated 
that for any usual source of matter the theory was not 
stable, but intrinsically dynamical p^ . Famously, he in- 
voked the cosmological constant - his "groflte Eselei" - 
to avoid the issue. This paper attempts to lay the concep- 
tual groundwork for an alternative first principles route 
to the physical interpretation of cosmological general rel- 
ativity, taking account of observational evidence that is 
immeasurably better now than it was in 1917. 

The plan of the paper is as follows. Preliminary defini- 
tions, motivations, and a statement of the cosmological 
equivalence principle are presented in Sec. |lTl In Sec. IIIII 
the thought experiments introduced in Sec.|TT]are gener- 
alised to the case of regions of different density, and an 



estimate of clock rate variance is given based on presently 
observed density contrasts. The definition of the cosmic 
rest frame, and its relation to other frames used in cos- 
mological averaging, is clarified in Sec. IIVI In Sec. |V] a 
numerical estimate of the time-varying relative decelera- 
tion of voids relative to walls, where galaxies are located, 
is computed over the lifetime of the universe, and its 
cosmological implications discussed. The role of initial 
conditions and a possible conceptual relationship to Pen- 
rose's Weyl curvature hypothesis are discussed in Sec.lVIl 
The paper concludes with a summary discussion, includ- 
ing some further speculations, in Sec. IVIll 



II. AVERAGING AND THE EQUIVALENCE 
PRINCIPLE 

The SEP addresses physics on relatively small scales 
in which matter can be treated by field theories in a lo- 
cal inertial frame (LIE). In a LIF fundamental nongrav- 
itational interactions are described by Lagrangian field 
theories, and in four spacetime dimensions the mass, m, 
and spin, s, are given by the Casimir invariants of rep- 
resentations of the local Lorentz group, P^P^ — —m?(? 
and W^W^ = s{s + l)m^c^ti^, where P^ is a particle 
4-momentum and the Pauli-Lubanski pseudovector. 
To make a transition to macroscopic scales we have to 
take averages of these quantities to obtain a hydrody- 
namic limit: a fluid description with an effective energy- 
momentum tensor, not derivable from a field Lagrangian. 

Following the end of the radiation dominated epoch, 
on cosmological scales of averaging, we largely deal with 
situations in which the massive matter fields form con- 
densates for which the average spin is zero, and can be 
treated as pressureless dust described by one scalar pa- 
rameter, the density, p. Neutrinos, on account of their 
ultralight masses, near relativistic speeds, and their in- 
teresting property that mass eigenstates do not coincide 
with flavour eigenstates, are an exception to this rule. 
However, neutrinos are so light that from the point of 
view of the cosmological background they can be effec- 
tively treated as a massless species. 

Massless particles which are relevant for cosmology will 
be considered. Gravitational waves will not be considered 
here, as their cosmological contribution at late epochs is 
negligible. Electromagnetic waves are considered, as they 
are the means by which almost all our information about 
the universe is transmitted. For the purpose of cosmo- 
logical averages it is sufficient that light propagation can 
be treated in the geometric optics limit, and that the cos- 
mic radiation background has a perfect fluid equation of 
state P — \pc^ . 

The question of the largest scale on which the SEP 
can apply has not, to the best of my knowledge, been ad- 
dressed in a fundamental way. Taken literally, as soon as 
we are dealing with scales on which particle masses must 
be treated as an integrated density, then the SEP can 
only apply if we ignore the gravitational interaction. Of 
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course, in practice, for many practical purposes we can 
treat the gravitational interaction by Newtonian gravity 
on the scale of the solar system, and the scale of galax- 
ies. However, the true nature of general relativity is to 
replace the Newtonian gravitational force by a dynam- 
ical spacetime, not to replace spacetime by Newtonian 
gravity in the limit of weak fields. 

It must be recalled that whereas Newtonian gravity 
is a nonlocal interaction on a static space, general rel- 
ativity is a local dynamical theory of spacetime. The 
curvature of the spacetime background is not local, and 
the Newtonian limit picks up the aspect of the nonlocal 
curvature of the background generated by slowly moving 
point particles, but in full general relativity, changes in 
the curvature of the background can only be built up over 
time by locally propagating processes. This means that 
even in the weak-field limit there exist circumstances in 
which Newtonian gravity on a static background is not 
an appropriate limit. In cosmology the background is 
not static, and thus clearly the Newtonian approximation 
may break down. While much has been achieved numer- 
ically by assuming TV-body Newtonian interactions on a 
pre-existing FLRW background, the universe is clearly 
inhomogeneous on scales of at least tens of megaparsecs 
at the present epoch. To consider the dynamical gravita- 
tional processes which lead to macroscopic variations of 
the spacetime background over these scales, we need to 
go back to first principles. 



A. The cosmological equivalence principle 



The rationale for the statement of the cosmological 
equivalence principle (CEP) is twofold. First, if we are 
to demand a smooth Newtonian gravitational limit in all 
circumstances, we have to deal with the fact that Newto- 
nian gravity deals with just one scalar source, the density, 
whereas general relativity is tensorial. This means that 
we must be dealing with an average spacetime with sym- 
metries in taking a Newtonian gravity limit. Since grav- 
itation with matter is dynamical, it is the symmetries 
involving time that should be removed from the SEP to 
account for the average density of matter. The metric (H]) 
does this while preserving the isotropy and homogeneity 
of space regionally within a cosmological inertial frame 
(GIF). 

Second, at the core of the equivalence principle is the 
notion that we can always choose a LIF, for example, by 
specifying Riemann normal coordinates. However, in the 
case of the volume-expanding and contracting motions 
of the metric ((!]), as illustrated by Fig. [1] it is impossible 
to locally distinguish between the case of comoving par- 
ticles at rest in an expanding metric ([1]) and the case of 
particles in motion in the static Minkowski space of the 
relevant LIF, a point which has been emphasised by a 
number of authors recently from various points of view 
On local scales, both yield the Hubble law redshift 
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My proposed answer to the question of the largest scale 
on which the SEP can apply is to deal with the aver- 
age effects of density by extending the SEP to poten- 
tially larger scales while removing the time-translation 
and boost symmetries of the background as follows: 

At any event, always and everywhere, it is possible to 
choose a suitably defined spacetime neighbourhood, the 
cosmological inertial frame, in which average motions 
(timelike and null) can be described by geodesies in a ge- 
ometry which is Minkowski up to some time- dependent 
conformal transformation, 

ds^^p = a^i-q) [-drf + dr^ + r'^d^f] , (1) 

where dfi^ is the metric on a 2-sphere. 

The standard SEP is obtained in the limit that a{rj) is 
constant, which physically corresponds to virialised sys- 
tems that can be effectively thought of as asymptotically 
flat. Alternatively, during very short time intervals over 
which the time variation of the scale factor can be ne- 
glected, the standard SEP is also retrieved. The idea 
here is that even when spacetime is dynamical in cos- 
mology, one can always find a spacetime neighbourhood 
whose average volume expansion can be characterised by 
a metric of the form ([T]) for time intervals over which airf) 
varies. The metric ([T]) is of course that of the spatially 
flat FLRW geometry in conformal coordinates. 



where ir is the radial proper distance from an observer at 
the origin to a source, and an overdot denotes a derivative 
with respect to t where cdt ~ adrj. This is true whether 
the exact relation, z + 1 = a^/a, is used or the radial 
Doppler formula, z + I = [{c + v)/{c — v)]^^^, of special 
relativity is used, before making a local approximation 
[T^ . Mathematically the equivalence of the two situa- 
tions might be viewed as a consequence of d/drj being a 
conformal Killing vector of ([T|). 

The aim of the CEP is to go beyond the limit of the 
static special relativistic LIF, to consider arbitrarily long 
time intervals over which the motion of the particles is 
decelerated, a < 0, by the average density of matter. 
As Einstein himself stated [l^l, "In a consistent theory 
of relativity there can be no inertia relatively to 'space', 
but only an inertia of masses relatively to one another" . 
Since the deceleration of the volume expansion is due 
to the backreaction of the average density of the matter 
particles in defining their average background, the CEP 
thus represents a refinement in the definition of inertial 
frames. To demonstrate this at a conceptual level, we will 
first show that for localised regions a suitable equivalent 
of decelerated Minkowski particles can always be found 
for the motion of a congruence of comoving particles in 
P]), even for arbitrarily long time intervals. 
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arbitrarily long supply of string is wound. The strings 
initially unreel at the same uniform rate, representing a 
"recession velocity" . Each observer carries synchronised 
clocks [l^, and at a prearranged local proper time all 
observers apply a braking mechanism, the braking mech- 
anisms having been pre-programmed to deliver the same 
impulse as a function of local time. 




FIG. 1: A set of particles undergoes an isotropic spatial 3- 
volume expansion in a spatially Bat local region. For the same 
initial conditions, provided that we consider time intervals 
over which deceleration of the expansion is negligible, it is 
impossible to locally distinguish the case of particles at rest 
in a dynamically expanding cosmological space from particles 
moving isotropically in a static Minkowski space. One spatial 
dimension is suppressed. 



FIG. 2: The semi-tethered lattice. Point particle observers in 
a homogeneous spatial lattice, initially expanding with uni- 
form velocity, are attached to nearest neighbours by strings. 
Each string is fixed at one end. The arrowheads denote the 
free end of each string wound on a spool at a neighbouring 
particle, to which the observers apply brakes in a synchro- 
nised fashion according to a pre-determined plan. The time 
evolution of the lattice follows a course similar to that of the 
spatial grid in Fig. [IJ but with deceleration. 



B. Thought experiment: The semi-tethered lattice 

Given a spatial point, taken to coincide with the origin 
of the coordinates ^ , we might be tempted to think that 
we can mimic the situation of an arbitrary decelerating 
expansion of ^ in Minkowski space by attaching rockets 
to a set of test particles initially chosen to be equidistant 
from the given centre and with identical radial velocities. 
The rockets could then be fired in unison at the same rate 
to give an inward radial acceleration while maintaining 
equidistance of the particles from the centre and equality 
of instantaneous radial velocities as seen from the frame 
of the central observer with global Minkowski time, t. 

The rocket analogy is flawed, however. To see the flaw, 
let us consider a second Minkowski space thought exper- 
iment which is directly equivalent to the cosmological 
situation of comoving particles moving in the geometry 
(III), even in the case of deceleration of the volume flow. 
Let us construct what I will call the semi-tethered lat- 
tice by the following means. Take a lattice of Minkowski 
observers, initially moving isotropically away from each 
nearest neighbour at uniform initial velocities [l^. The 
lattice observers are chosen to be equidistant along mu- 
tual oriented i, y and z axes. Now suppose that the 
observers are each attached to six others by strings of 
negligible mass and identical tension along the mutually 
oriented spatial axes, as in Fig. [21 The string in each 
observer's negative x, y and z directions is held fixed 
and extends to the observer's nearest neighbour in those 
directions. The string extending towards each nearest 
neighbour in the positive x, y and z directions unreels 
freely from a spool at the observer's site on which an 



The semi-tethered lattice is directly analogous to the 
case of the decelerating volume expansion of ((T]) due to 
some average homogeneous matter density, because it 
maintains the homogeneity and isotropy of space over 
a region as large as the lattice. In the case of the fir- 
ing of rockets, the act of firing a rocket means that each 
observer with a rocket feels a net force in a particular 
direction, while also breaking the symmetry of the homo- 
geneity of space. In the case of the semi-tethered lattice 
work is done in applying the brakes, and energy can be 
extracted from this. However, since brakes are applied 
in unison (according to local proper time at each lattice 
site), there is no net force on any observer in the lat- 
tice. Although the rate at which the brakes are applied 
can be an arbitrary pre-arranged function of local proper 
time, provided the braking function is applied uniformly 
at every lattice site, the clocks will remain synchronous 
in the comoving sense, as all observers have undergone 
the same relative deceleration. 

The semi-tethered lattice is also a useful analogy be- 
cause the kinetic energy of the particles is converted di- 
rectly to heat in the brakes, which might then be con- 
verted to other forms. It is a direct analogue for the 
conversion of kinetic energy of particles in the expansion 
of the universe to other forms of energy via gravitational 
collapse. Apart from the energy of massless and near 
massless species released during recombination and ear- 
lier phase transitions, all the useful forms of energy in 
the present observable universe have undergone such a 
process of transformation: from the kinetic energy of ex- 
pansion to gravitational energy, and then to other forms. 

The semi-tethered lattice and the CEP might be said 
to extend the elevator thought experiments and the 1907 
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Einstein equivalence principle in a natural fashion to 
the case of a homogeneous isotropic nonstatic gravita- 
tional field [l^. The important point in the present 
situation is that in both the cosmological case and the 
semi-tethcrcd lattice analogy the observers feel no net 
force from the relative deceleration, which justifies the 
description of ([T]) as a cosmological inertial frame. 



C. Cosmological inertial frames and averaging 

Although ll]) is simply the standard spatially flat 
FLRW geometry, the important physical difference from 
the usual treatment of cosmology is that here it is not to 
be viewed as a global metric of spacetime, but as an av- 
erage cosmological frame over some region. The FLRW 
geometries with spatial sections of positive and negative 
spatial Gaussian curvature, k, 



FLRW 



-d^ 



(2) 



can of course also be brought into the form ([T]) over small 
spatial domains as compared to the radius of curvature. 
This follows directly from the fact that the FLRW geome- 
tries are conformally flat and coordinates can be found 
to bring ([2]) into the form 



ds^ 



^ a^{r],r)[-dr]^ + dr^ +r^dn^] . (3) 



The new physical interpretation is that no single 
FLRW geometry is to be taken as a global average 
geometry for the whole universe for all times. However, 
with an extension of the SEP to account for the cosmo- 
logical average effects of the density of matter, regional 
cosmological inertial frames with average geometry ([T]) 
can always be found. It should also be emphasised that 
since these are average frames, with differing regional 
scale factors and local coordinates, no metric of the form 
(12) is substituted into the Einstein equations and solved. 
Rather, an appropriate average of the full inhomogeneous 
Einstein equations, such as a Buchert average [17], should 
be applied to solve for the background average of the 
inhomogeneous geometry. One can solve the Buchert 
equations for a realistic approximation to the observed 
universe 0, but care must be taken in interpreting the 
solution, as we must account for where the observers are 
within the inhomogeneous structure when it comes to the 
relative calibration of their rods and clocks [ij . 

In the case of a semi-tethered lattice which is confined 
to a finite region, the relative deceleration of the region 
would give a local proper time to comoving observers 
in the lattice different from that of a global Minkowski 
observer, even if this observer had a synchronous clock 
before the volume deceleration began. By the CEP the 
average homogeneous density in different regions likewise 
sets a standard of local time for CIFs, and this may vary 
significantly when regional density contrasts grow. 



Since any GIF ([T]) is an average, the relationship of 
the average to the inhomogeneous geometry needs to 
be stated. For example, in terms of the proper time, 
cdt = adrj, of an observer "comoving" with the GIF, the 
uniform expansion of the GIF should be viewed as an 
average only: 



GIF 



H. 



GIF' 



(4) 



where 6 is the volume expansion, angle brackets de- 
note the appropriate average, and an overdot denotes a 
derivative with respect to t. In relation to the averaging 
scheme, the specification of a GIF should incorporate a 
notion of "centre-of-mass motion" in the sense that the 
variance in the GIF volume expansion is an average of 
shear and vorticity fluctuations for which the net backre- 
action within a GIF vanishes. In the notation of Buchert 
and Garfora llSl. 



GIF 



GIF 



(^)gif 



GIF' 



(5) 



in Buchert 's scheme with vanishing vorticity, where — 
^cTfiu'y'^'^ is the scalar shear. The condition (O ensures 
that the average kinematic backreaction within a GIF 
vanishes in Buchert 's scheme. For GIFs containing galax- 
ies, vorticity will be important, and so Buchert's scheme 
needs to be generalised to include vorticity before a pre- 
cise statement can be made. As the contributions of vor- 
ticity and shear are of opposite sign in the Raychaudhuri 
equation their average contributions may even be largely 
self-cancelling, for virialised systems at least. However, a 
precise formalism including vorticity remains to be devel- 
oped. It is possible that ideas could be incorporated from 
Zalaletdinov's averaging scheme [l^, though the physi- 
cal content of that scheme first needs to be more clearly 
elucidated in a cosmological setting [2^ . 

For the present paper it is sufficient that the back- 
reaction can be neglected in characterising the average 
properties within a GIF. Thus a GIF is characterised by 
a single scalar, the volume expansion In an expand- 
ing universe with collapsing substructure, this volume 
expansion is best viewed as a macroscopic property of a 
given GIF, but often not of its more finely grained sub- 
regions. In particular, if a GIF contains a galaxy or a 
galaxy cluster it will contain virialised regions and may 
also contain collapsing regions. For cosmological aver- 
ages, when one is interested in comparing deceleration 
rates in expanding regions of different density, the no- 
tion of a finite infinity cutoff scale [U, [2l[ is suggested as 
a minimum region for a GIF in relevant averages. 



III. THOUGHT EXPERIMENT: RELATIVE 
HOMOGENEOUS ISOTROPIC DECELERATION 

The beauty of the equivalence principle is that it allows 
one to quantitatively deduce the order of magnitude of 
simple effects, such as the leading order of gravitational 
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redshift [1], by simple thought experiments alone. The 
basic physical effect that is of interest here - the gravi- 
tational energy cost of a spatial curvature gradient - can 
likewise be understood by a simple thought experiment. 

From the evidence of the cosmic microwave back- 
ground (CMB) radiation, we know that, apart from tiny 
fluctuations of order 5plp^lO~^ in photons and the 
baryons they couple to, and density fluctuations perhaps 
an order of magnitude larger in nonbaryonic dark mat- 
ter particles, the observable universe was close to spa- 
tially flat, homogeneous, and isotropic at the epoch of 
last scattering. Assuming the Copernican principle, it 
was sufficient to describe the universe by a single frame 
(P) at that epoch. However, since a GIF is a local re- 
gional frame, we should be careful never to construct a 
GIF with spatial extent larger than the particle horizon 
at any epoch. Rather, it is better physically to think of 
a class of local GIFs at disjoint spatial locations at the 
surface of last scattering, in which there exist geodesic 
congruences of observers like those depicted in Fig. [U 
which all have an identical uniform expansion away from 
each other, on account of the initial uniformity of the 
Hubble flow at that epoch. 

Let us first analogously consider two sets of disjoint 
semi-tethered lattices, with identical initial local expan- 
sion velocities, in a static background Minkowski space. 
[See Fig. [3)^a).] The observers in the first congruence 
apply brakes in unison to decelerate homogeneously and 
isotropically, with inward 4-acceleration of magnitude, 
a^{T^), as measured by the force applied to the brakes in 
the frame of any of the observers, where is the proper 
time measured by any of them. From the viewpoint of 
a global Minkowski observer, members of the congru- 
ence will agree on their measurements of time, t^^, though 
this time of course differs from the global Minkowski ob- 
server's time, t. 

Now take a second semi-tethered lattice, with the same 
initial expansion speed, where brakes are applied with 
a force corresponding to a 4-acceleration of magnitude 
a^{T^). At any global Minkowski time, t, we will as- 
sume that when transformed from their proper frames to 
that of the global Minkowski observer, at each time step 
a^{t) > a^{t) > 0. It is then the case that the members 
of the first congruence decelerate more than the mem- 
bers of the second congruence, and at any time t the 
proper times satisfy <- '''2- The members of the first 
congruence age less quickly than members of the second 
congruence. 

By the CEP, the case of volume expansion of two re- 
gions of different average density at late times is entirely 
analogous. The fact that we are able to apply the equiv- 
alence of the two circumstances rests on the fact that 
the expansion of the universe was extremely uniform at 
the time of last scattering, by the evidence of the GMB. 
At this epoch all regions effectively have the same den- 
sity - apart from negligible fluctuations ~ and the same 
uniform Hubble flow. At late epochs, suppose that in 
the frame of any average cosmological observer there are 



less deceleration 



more deceleration 



(a) 



average t = const 



^ ^ more dense 



less dense 



gradient in <R> 



(b) 



last-scattering 



FIG. 3: Two equivalent situations: (a) in Minkowski space 
observers in separate semi-tethered lattices, initially expand- 
ing at the same rate, apply brakes homogeneously and isotrop- 
ically within their respective regions but at different rates; 
(b) in the universe which is close to homogeneous and 
isotropic at last-scattering, observers in separated regions ini- 
tially recede from each other at the same rate, but experience 
different locally homogeneous isotropic decelerations as local 
density contrasts grow, fn both cases there is a relative decel- 
eration of the observer congruences, and those in the region 
which has decelerated more will age less. 



regions of different density which have decelerated by dif- 
ferent amounts by a given time, t, according to that ob- 
server. Then, by the CEP the local proper time of the 
isotropic observers in the denser region, which has de- 
celerated more, will be less than that of the equivalent 
observers in the less dense region which has decelerated 
less. [See Fig. [21b).] Consequently the proper time of the 
observers in the more dense CIF will he less than that 
of those in the less dense CIF, by equivalence of the two 
situations. 

The fact that a global Minkowski observer does not ex- 
ist in the second case does not invalidate the argument. 
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The global Minkowski time is just a coordinate label, 
and in the cosmological case the only restriction on the 
t =const slices is that the expansion of both average con- 
gruences must remain homogeneous and isotropic in local 
regions of different average density in the global t =const 
slice. Of course, we need to be careful to patch together 
different CIFs continuously to specify the slice, as we 
will further discuss in Sec. IIVI In this way the equiva- 
lence to the Minkowski space case is maintained. Thus 
in the cosmological case, provided that we refer to local 
homogeneous isotropic expansion in different regions on 
any average t =const slice, (where t is some coordinate 
label), then if such regions are still expanding and have 
a significant density contrast, we can expect a significant 
clock rate variance. 

This equivalence directly establishes the idea of a gravi- 
tational energy cost for a spatial curvature gradient^ since 
the existence of expanding regions of different density 
within an average t =const slice implies a gradient in the 
average Ricci scalar curvature, (TZ), on one hand, while 
the fact that the local proper time varies on account of 
the relative deceleration implies a gradient in gravita- 
tional energy on the other. 



A. Order of magnitude estimate of clock rate 
variance 

An order of magnitude estimate of present epoch clock 
rate variances due to gravitational energy gradients in- 
duced by relative volume deceleration of the background 
can now be made by directly using observationally mea- 
sured density contrasts. Although a GIF is a frame ([1]) 
within which backreaction can be neglected, to determine 
its scale factor over long periods of time one must con- 
sider the evolution of the universe within which the GIF 
is embedded. Such evolution will, in general, include the 
effects of backreaction. However, if the backreaction is 
small, an order of magnitude estimate of the clock rate 
variance can be made assuming that regions with ob- 
served strong density contrasts evolve independently by 
solutions of the local Friedmann equation for regions of 
different density. There will be a relative deceleration of 
the local background of such regions, which via equiva- 
lence to the Minkowski space semi-tethered lattices, will 
accumulate clock rate differences. 

Galaxies formed from perturbations which were greater 
than critical density and if space is negatively curved on 
average, they must always be bounded by a region which 
is spatially flat. These on-average spatially flat locally 
expandirig bounding regions are called finite infinity re- 
gions P, im , and a union of such regions is called a wall. 
Since they are spatially flat, neglecting backreaction, the 
evolution of the wall GIFs can be approximated by spa- 
tially flat Einstein-de Sitter regions with local scale fac- 



tor a„= a, 



where Hi is the common initial 



Hubble parameter at last scattering, and is a constant. 
On the other, hand GIFs within voids can be approx- 



imated by portions of spatially open FLRW solutions, 
given parametrically by 



'^'^^ (coshTj — 1) , 



2(1 - n,) 



(sinh r] — rj) , 



(6) 



(7) 



where cdr^ = Ovdrj, fli — 1 — is an initial density 
parameter at last scattering, 1 and Oi is a constant. 
Using ^ and ([7]) in the Friedmann equation, one obtains 
a standard parametric relation for the density parameter. 



2 (cosh 77 — 1) 
sinh'^ rj 



(8) 



which is to be viewed here as a regional parameter in a 
GIF inside a void. _ 

We now follow the analysis of Ref . [22] , where the au- 
thor's proposal was first advanced. There an attempt was 
made to estimate cosmological parameters by making the 
approximation that the evolution of the observable uni- 
verse was entirely due to the voids. In fact, backreac- 
tion between the walls and voids must be included to 
obtain more reliable estimates of cosmological parame- 
ters [U, d^. However, if we confine attention to small 
regional GIFs, then the argument of Ref. [22] can give 
an estimate of clock rate variance from observed density 
contrasts. In particular, since the critical density also de- 
fines the Einstein-de Sitter standard of time of the wall 
GIFs, it also follows that 



o-vj ri?(coshry 



IP 



Gombining ([5]) and ^ we find 



fli{coahri — 1)^ 
r^i)'^/^ sinh rj 



3(1 



(9) 



(10) 



Differentiating both ^ and (fTO|) we find a relative clock 
rate, which we will call the lapse function, given by 



7(^ 



dTv 3(cosh77 + l) 
dTw 2(coshr7 + 2) 



(11) 



A relative clock rate variance due to the relative vol- 
ume deceleration between GIFsJn walls^and voids can 
now be estimated since ft — l-\-S, where S is the density 
contrast, and inverting ([8]) we have 



cosh 7] 



1 



1 



(12) 



Hoyle and Vogeley estimate that 40%~50% of the 
present-epoch universe is in voids of diameter about 
30h~^Mpc, with the statistics summarised in Table HI 
The density contrasts quoted are an average, and are of 
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greater magnitude in the centres of the voids which are 
extremely empty. Taking these values as indicative, if we 
assume S = -0.9 then ([TT]) and ([12]) give 7 = 1.42, while 
if we assume 6 = —0.95 then 7 = 1.46. For larger density 
contrasts, in the centre of the void the lapse would ap- 
proach the limiting value 7 1.5, which represents the 
relative local expansion rates of an empty Milne universe 
to an Einstein-de Sitter one. 
These clock rate variances of 42%-46% are large, and 
counter intuitive given we usually encounter large clock 
rate differences only for large local boosts or for den- 
sity contrasts from extremely compact sources in static 
backgrounds, such as black holes. However, the effect 
described is neither of these familiar situations, as the 
background is not static. Rather, it is the clock rate vari- 
ance due to the cumulative effect of a very small relative 
deceleration of the background. The above variances are 
simply those demanded by the CEP taking present epoch 
density contrasts observed in the actual universe. 



Survey 


Void diameter 


Density contrast 


PSCz 


(29.8 ± 3.5)/l- 


-^Mpc 


(5 = -0.92 ±0.03 


uzc 


(29.2 ± 2.7)/l" 


'^Mpc 


5 = -0.96 ±0.01 


2dF NGP 


(29.8 ± 5.3)/i" 


'^Mpc 


5 = -0.94 ±0.02 


2dF SGP 


(31.2±5.3)/i" 




5 = -0.94 ±0.02 



TABLE I: Dominant void statistics in the Point Source Cata- 
logue Survey (PSCz), the Updated Zwicky Catalogue (UZC), 
and the 2 degree Field Survey (2dF) North Galactic Pole 
(NGP) and South Galactic Pole (SGP), from Refs. [H. 

Taking backreaction into account to determine a full 
present epoch average [l|, 0], the present epoch lapse 
difference between walls and a volume-average loca- 
tion in a void based on luminosity distance data fits is 
7 = l.SSto^os @- This is an average value; the rela- 
tive lapse would be larger in the void centres. Thus the 
estimates made without backreaction are reasonably ac- 
curate, showing that the effect is not a direct consequence 
of backreaction in the evolution equations but rather of 
relative volume deceleration alone. From (|11|1 and p2|l 
without backreaction the density contrast estimate for a 
local GIF at the volume-average position is 5 = —0.83 
for a lapse of 7 1.38. If such a clock variance were pro- 
duced by a uniform acceleration in Minkowski space, a 
simple calculation shows it would require a tiny relative 
acceleration of order 5.5 x 10~^°ms~^ over the lifetime 
of the universe. Of course, such a relative acceleration is 
not uniform: a better estimate is presented in Sec. fVl 

The argument above relies on it being possible to 
choose a locally uniform Hubble flow gauge, as will be 
discussed in the next section. Such a gauge can be main- 
tained outside finite infinity regions, but not within them 
where collapsing regions are located, and where vorticity 
and tidal torques become important. Thus there are no 
obvious inferences analogous to those above that can be 
made with respect to bound structures from observations 



of the magnitude of their positive density contrasts. 

IV. AVERAGE FRAMES 
A. The cosmic rest frame 

In taking cosmological averages with inhomogeneous 
structures, the question arises as to which average frames 
have the most utility. One must generally make a choice 
of gauge in specifying such frames. One way of viewing 
the SEP is that we can always set the first derivatives 
of the metric to zero near a point. In particular, the 
volume expansion 9, which involves first derivatives of 
the metric, can always be set to zero. The CEP extends 
this by the statement that in the dynamical situation a 
gauge can be chosen in which the volume expansion in a 
CIF is spatially uniform, but varying with time. 

As was pointed out in Sec. lIIIl in order to compare CIFs 
in regions of different densities, we need to specify suit- 
able spacelike slices '26^ by patching together different 
CIFs in a continuous manner. Operationally, the way to 
do this is first to choose an orientation of the 4-velocity 
fields, d/dTj, of comoving observers in a CIF such that 
the CMB radiation is isotropic in each frame. In terms 
of the local proper time, t^, of such observers the metric 
^ is rewritten 

where the index / runs over CIFs. Secondly the spacelike 
slice is specified by the demand that the locally measured 
value of the volume expansion remains uniform as one 
moves from the patch of one CIF to the next. In other 
words, the "local Hubble fiow" remains uniform in this 
gauge even though the proper lengths and proper time 
scales will change as one moves between CIFs of different 
density. As discussed in Ref. [H, although the proper 
volume of voids increases faster than that of wall regions, 
this is compensated for locally by the faster relative rate 
of the void clocks. Relative to any one set of clocks, such 
as our own, it will always appear that voids expand faster 
than walls. So the average Hubble flow over both walls 
and voids - by one set of clocks - will generally differ 
from the underlying uniform flow. Its value is a choice of 
gauge depending on the choice of fiducial observer. 

The uniform Hubble flow slices defined in this man- 
ner constitute the cosmic rest frame: surfaces within 
which the CMB is isotropic, even though the mean value 
of the CMB temperature, and angular scale of CMB 
anisotropics, will vary from point to point as spatial vol- 
umes vary in relation to proper radius with changes in 
spatial curvature. The proper times of CIF observers 
within the slice will also vary. We can choose the clocks 
of a canonical set of observers in expanding regions with 
the same average local density to label the slices, pro- 
vided we realise that this time labelling is only a proper 
time in particular locations on the slice and is just a co- 
ordinate label elsewhere. 
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The uniform Hubble flow gauge is one of the standard 
gauges of perturbation theory in FLRW models ^] , and 
has been further refined with the addition of a minimal 
shift distortion condition by Bicak, Katz and Lynden- 
Bell [9]. These authors recognise the resulting "Mach 
1 gauge" as one of three possible gauges which best in- 
corporates Mach's principle, and within which there is a 
minimal amount of residual gauge freedom. Bicak, Katz 
and Lynden-Bell work within the framework of pertur- 
bations of a global FLRW geometry. The viewpoint of 
the present paper is that there is no single global FLRW 
geometry. However, on account of the CEP the spatially 
flat and negatively curved FLRW geometries can be con- 
sidered as local regional geometries on spacelike slices as 
the density varies relative to the critical density. Thus a 
modiflcation of the formalism of Bicak, Katz and Lynden- 
Bell may be an appropriate starting point to deal with an 
averaging formalism for the full nonlinear inhomogencous 
problem, as an alternative to Buchert's scheme. 



B. Buchert averaging 



Buchert's averaging scheme [17[ is based on the start- 
ing point that, in the case of an energy-momentum tensor 
for irrotational dust particles in the presence of inhomo- 
geneities, one can choose Gaussian normal coordinates 



(14) 



comoving with the dust. The scalar density appearing in 
the energy-momentum tensor, 



where n'^ 



dt 



, then represents the rest mass density of 
the dust, and one averages over spatial slices of constant 
t orthogonal to the flow, over regions which conserve the 
rest mass of a portion of the fluid in a domain, P, with 
continuity equation 



a 



dt{p) + ^{p) = Q 



a 

,1/3. 



(15) 



where a(i) = \V{t) /V{t^)\' vfit\iV{t) = ^^d^x^/det^. 
Here angle brackets denote the spatial volume average of 

a quantity, so that (7^) = d^xy/det TZ{t, x)) /V{t) 

is the average spatial curvature, for example. The 
Buchert equations consist of and 



3a 



^ = &ttG{p) - ic2(7^) - iQ, 



- = -47rG(p) + Q, 
a 

9t(a'Q)+aVat(a2(7^)) =0, 



(16) 

(17) 
(18) 



where Q = ^ {{0'^) - (6)'^) - 2{a'^), is the kinematic 
backreaction. Equation psp is an integrability condition 
which ensures closure of the other equations. 



Since the backreaction term Q includes variance in vol- 
ume expansion, and this is to be evaluated on a constant t 
slice, it is clear that as compared to the cosmic rest frame 
of Sec. IIV Al different physical premises underlie the in- 
terpretation implicitly assumed by Buchert's scheme. My 
approach is therefore different from other approaches to 
cosmological building that have been adopted in the con- 
text of Buchert averaging 28, 29, 30, 31]. The differences 
may be understood from the fact that in talking about 
the "rest mass density of the dust" one is actually dealing 
with a concept which depends on the manner in which 
dust "particles" are coarse grained. Since all forms of 
energy have a rest mass equivalent, the kinetic energy of 
particles within a dust particle is included as rest mass. 
Similarly, since Ricci curvature affects spatial volumes 
relative to their diameter, the concept of a rest density 
depends on the scale of coarse graining relative to the 
curvature scale. 

In general, the notion of "comoving with the dust" 
implicit in Buchert's scheme can be very distinct from 
"comoving with the background" , although the notions 
coincide for FLRW models. This is well illustrated by 
the exact spherically symmetric Lemaitre-Tolman-Bondi 
(LTB) models [s^l for pressureless dust, with a prescribed 
inhomogencous density p(t,r). These can be written in 
Gaussian normal form (I14p , making it straightforward to 
compute a Buchert average |3l, [l^l [111 . At fixed comov- 
ing proper time, t, as the radial coordinate r varies the 
LTB dust shells have different densities, different spatial 
curvature, nonzero shear, and, in general, observers at 
r > would not expect to see an isotropic CMB. Since 
the solution is completely dynamical, there is no average 
homogeneous isotropic background with respect to which 
one could be comoving, unless one puts in such a back- 
ground by hand by making the model asymptotic to an 
FLRW model at large r. 

With respect to fixed FLRW backgrounds, an alterna- 
tive simple way to treat spherical inhomogeneities is by 
the spherical top hat model, using concentric spherical 
shells [36|, |37|. In the case of a void in a background 
Einstein-de Sitter universe, for example, a spherical un- 
derdense shell will acquire a peculiar velocity with respect 
to the background which tends to 50% of the background 
Hubble rate at late times One can account for the 
kinetic energy of the shell, but in view of the large pecu- 
liar velocity of the shell there is a limit to the extent to 
which it can be considered comoving with respect to the 
background with a synchronous clock. 

Once one averages on the scale of statistical homogene- 
ity, as in (|15p . one wants to have a sense of "comoving 
with the background"; i.e., different observers in differ- 
ent averaging cells should have a notion of determining 
the same average density at the same cosmological epoch, 
and one should be able to talk about motion with respect 
to canonically defined observers. 

In general, when the background is only statistically 
homogeneous and isotropic, there is an ambiguity in dis- 
tinguishing between motion of the background and mo- 
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tion with respect to the background. The spirit of the 
CEP is that in the case of volume expansion, this is be- 
cause there is a fundamental indistinguishability: the 
Hubble parameter is a gauge choice. The definition of 
the uniform Hubble expansion gauge of Sec. lIVAl makes 
an unambiguous separation of the "kinetic energy of the 
volume expansion of space" for regions which are locally 
spatially flat from other forms of energy. In the stan- 
dard interpretation of Buchert averaging, depending on 
the choice of averaging volume and the manner in which 
one coarse grains over dust cells, the kinetic energy of 
the volume expansion of space can be intermingled with 
other forms of energy. 

The view that I adopt is that one can either use 
Buchert averaging in defining CIFs on small scales by the 
requirement that the kinematic backreaction can be ne- 
glected, as in (O, or alternatively, on dust cells of at least 
the scale of statistical homogeneity, of about 100/i~^Mpc, 
at which scale the volume expansion is statistically uni- 
form. In both cases the time parameter and {TZ) can 
be regarded as relevant parameters which describe the 
collective degrees of freedom of the cell. However, they 
should not be regarded as representative for observers 
in finer partitions within the cell. For CIFs containing 
galaxies with black holes this imperative is obvious. In 
the case of statistically homogeneous cells in cosmology, 
the case has been less obvious, but I believe it is equally 
imperative on account of the density contrasts that are 
observed between voids and walls below the scale of ho- 
mogeneity, together with the arguments of Sec. IIIII 

I will take the view that the Buchert time parameter is 
the relevant one for an observer at a volume-average po- 
sition in a statistically homogeneous cell. As the present 
universe is dominated by voids, this will be in a void 
but not at the void centres. Kinematic backreaction be- 
tween voids and walls from the volume-average perspec- 
tive must be included in determining the average evolu- 
tion. Einstein's equations are causal and depend on all 
events within the past light cone. Thus some sort of spa- 
tial averaging, such as Buchert averaging, is required to 
determine cosmic evolution. Buchert 's equations should 
thus be viewed as evolution equations. 

Observations are made on null cones, however. Thus a 
Buchert average is not the one we perform operationally 
in determining cosmic averages. Instead, a radial null 
geodesic average of a solution to the Buchert equations, 
combined with an operational identification of relevant 
classes of observers within the inhomogeneous structure, 
is required to make comparisons with observations. This 
is the approach adopted in Refs. [H, 



V. RELATIVE DECELERATION OF THE 
BACKGROUND AND THE WEAK-FIELD LIMIT 

In Refs. [J Q a model universe was constructed based 
on a regional division of cells of average homogeneity 
into voids, and the bubble walls that surround them, as- 



suming that backreaction within the walls and voids can 
be neglected, but not in the combined average [sol. [ioj. 
Space within the walls is assumed to be spatially flat on 
average, and space within the voids is negatively curved. 
Technically a "wall" is understood to be a union of con- 
nected finite infinity regions Q, namely CIFs in regions of 
average critical density. Observationally, the walls would 
include all morphological types of extended structures 
containing galaxy clusters: namely "sheets", "filaments" 
and "knots" [4l|. 

Qualitatively the regional division into voids and walls 
may well be a consequence of the evolution of the scales 
with statistical excesses of power in the primordial angu- 
lar power spectrum. The lOOh^^ Mpc scale of statistical 
homogeneity would correspond to the first Doppler peak 
(BAO scale) : on account of the finite sound speed in the 
primordial plasma no structures in excess of this scale 
are expected statistically, with the exception of those 
random structures that arise by percolation. The scale 
of the 30h~^ Mpc dominant voids would correspond to 
the evolution of the second Doppler peak, namely, the 
first rarefaction peak, which is well within the nonlinear 
regime of structure formation. The third Doppler peak, 
which is the first compression peak within the nonlinear 
regime, would give the scale of the largest bound struc- 
tures that have broken from the Hubble flow, namely, 
clusters of galaxies. Finite inflnity represents a demar- 
cation scale of on-average spatially flat regions between 
clusters of galaxies and voids. The fourth Doppler peak, 
the next rarefaction peak, may possibly give an indepen- 
dent scale corresponding to mini voids. In the two-scale 
approximation [42] of Refs. [U, 0|; we assume that the av- 
erage curvature of minivoids and dominant voids is the 
same. A further refinement might separate these scales. 
Ultimately these qualitative speculations about the cor- 
respondence of the Doppler peaks to the observed scales 
of present epoch structures should be verified from a nu- 
merical model of structure formation. 

An exact solution to the Buchert evolution equations of 
the two scales was found 01 , which by a Bayesian analy- 
sis fits the RiessOT gold supernovae data set [ij] at a level 
which is statistically indistinguishable from the standard 
spatially flat ACDM model _23||. The same best-fit pa- 
rameters also fit the angular scale of the sound horizon 
seen in CMB data, and the effective comoving baryon 
acoustic oscillation scale seen in angular diameter tests 
of galaxy clustering statistics [ll, [1^ . 

Those sceptical of the proposal sometimes question 
how the relatively large present lapse of 7 = 1.38 be- 
tween wall observers and the volume average can have 
possibly arisen. As pointed out already in Sec. IIII Al a 
relatively small relative deceleration of the background 
for the lifetime of the universe is sufficient to achieve 
this. Since the accumulated average lapse function 7 is 
not uniform in time, the equivalent relative deceleration 
of the background is not uniform. In this section we will 
estimate the relative deceleration which would produce 
the tracker solution mean lapse function 0]. 
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In principle, one is trying to compare the relative de- 
celeration of regions of different density as a function 
of the expansion age of the universe of some fiducial ob- 
server such as ourselves. Ideally, if sources did exist freely 
falling in voids unbound to condensed structures, then by 
the assumptions of Ref. [H they would have different red- 
shifts to sources in bound structures that coexist "at the 
same epoch" on account of the accumulated gravitational 
energy differences. We would have 



1 + 



(-k-Uw) 
(-k-Uob.) 



(19) 



for the redshift of a wall source with 4-velocity as 
seen by an observer with 4- velocity Uobs i where k is the 
4-velocity of a radial null geodesic. Similarly 



1 



(-k-U.) 
(-k-Uobs) 



(20) 



where Uv is^ a volume-average void source, and ^ 
In general [ij] 



1 -I- Zw = (l-f 2v)/7 



(21) 



where similarly to ([TT|) di = 7 drw represents the relative 
clock rates of volume-average observers to wall observers 
at the epoch of emission. Unfortunately, we have a mass 
biased view of the universe and only observe sources in 
bound systems in regions which are locally spatially flat 
on average; otherwise the relative blueshift —1 -I- 7"^ of 
a volume-average void clock relative to a wall one would 
have an obvious observational signature nearby. 

To determine the relative deceleration, we will invoke 
the CEP by considering the Minkowski space equivalent 
semi-tethered lattice analogy of Fig. [31 Ideally, we should 
have to calculate the difference in the rate of extraction of 
energy in applying brakes at different rates, to represent 
regions of different density in the actual universe. How- 
ever, this should be equivalent to asking what relative 
volume deceleration would be required to produce 7(tw) 
at any epoch, if 7(tw) is treated as a Lorentz gamma- 
factor in special relativity, beginning from two regions 
with the same initial expansion velocity (which is the 
situation at last scattering). In relation to (HU), this is 
equivalent to treating the blueshift —1 -1-7"^ of voids rel- 
ative to walls as if it were a standard transverse Doppler 
shift in special relativity. Since we are dealing with an 
isotropic volume deceleration there is no directional sig- 
nificance associated with a "direction of motion" in the 
special relativistic analogy. 

We assume equivalence to the special relativistic 4- 
acceleration a = where f/'' = 7(0,'!;*), which has a 
magnitude Q 
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with dt — ^ dr. Of course, in the present case we are 
not really dealing with a proper acceleration in a sin- 
gle direction, as the appropriate analogy is that of two 



semi-tethered lattices in which all directions contribute. 
Assuming equivalence of the situations the relative accel- 
eration of the background here has a magnitude 
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where, following the notation of Refs. P, 0]j i is the time 
parameter of the volume-average observer in a void, and 
dt = 7dTw, where is the time for an observer at finite 
infinity in a wall, which will be close to the time in an 
average galaxy. Furthermore, H is the bare or locally 
measured Hubble parameter, while on account of eq. (42) 
of Ref. %H = ^H — 7, is the dressed Hubble parameter 
obtained by averaging over both walls and voids on the 
past light cone. 

Using the tracker solution of Ref. 0] , 
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where is the present epoch value of the bare Hubble 
parameter B{t), fvit) is the void fraction and /vo its 
present epoch value. The void fraction here is that of all 
voids, including mini voids [45'| , and not just the dominant 
voids of Table HI The dressed Hubble parameter satisfies 
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while the time parameter t„ of wall observers is related 
to that of volume-average ones by 
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where f2 
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(1 ^ /vo)(2 -I- /vo) is the present epoch 



dressed matter density. From ([23l) -(p6 |l we find 
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Estimate of relative deceleration scale 



Using the estimates of Ref. j23l l. at the present epoch 
the void fraction is /vo = 0.76to'o9: while the dressed 

7+1-2 1 „„„-liv^ -1 



Hubble constant is = 61.7li['^km sec~^Mpc~ 
where the uncertainties are la uncertainties from a fit to 
the Riess07 gold data set The bare Hubble constant 
is then = 48.2^2 4 km sec"^Mpc~-^. From these val- 
ues the present epoch magnitude of the relative decelera- 
tion (|29p is ttp = 6.7l34xlO~^^ms~^. Furthermore, since 
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FIG. 4: The magnitude of the relative deceleration scale, a, 
as a function of redshift: (a) for redshifts z < 0.25; (b) for 
redshifts z < 2. The central curve shows the value for the 
best-fit parameters /vo ~ 0.76, = 61.7 km sec~^ Mpc~^ 
(Hg = 48.2 km sec~^ Mpc'^) from Ref HI/. The dashed 
curves show the corresponding results with la uncertainties, 
which are largely due to the large uncertainty on /vo, which 
is not tightly constrained by supernovae data. The upper 
dashed curve corresponds to /vo ~ 0.67 and the lower dashed 
curve to /vo = 0.88. In panel (a) the horizontal dotted 
lines indicate the bounds of the empirical acceleration scale of 
MOND when normalised for H^^ = 61.7 km sec~^ Mpc~^ . 



q; is a time-varying quantity, its best-fit value is plotted 
as a function of redshift for recent epochs in Fig. 21 

It is interesting to note that over the range of redshifts 
z < 0.25 in Fig. IH the curve for the best-fit parameters 
with present epoch void fraction /vo = 0.76 precisely tra- 
verses a range of values for a that has been used for the 
empirical acceleration scale of the modified Newtonian 
dynamics (MOND) scenario l4q|. In particular, using a 
range of recently quoted values [43] , the MOND acceler- 
ation scale is Ofmond — 1-21q'2 X I0~^^hjgms~^ , where 
/i75 = _ff(,/(75km sec~^Mpc~^). Using our best fit 
iJp = 61.7li['^km sec^^Mpc^^, this range corresponds 
to a„,o„d = S.ltl l X 10-"ms-2. 

The fact that a is larger at higher redshifts reflects the 



property that it scales in proportion to the bare Hubble 
parameter, H{t), as can be seen from (f^ . Both the 
bare Hubble parameter and dressed Hubble parameter 
are of course larger at earlier epochs at higher redshift, 
since the universe is always decelerating in the present 
model. (As demonstrated in Refs. [11, Q cosmic accelera- 
tion is a purely apparent effect related to clock rate vari- 
ance.) From one sees that a also scales in proportion 
to fv{t) which is smaller at earlier times, contributing a 
term in competition with both the term H(t) in the nu- 
merator of (|29p and with terms in the denominator. As 
a proportion of the Hubble flow at any epoch, the rela- 
tive deceleration is suppressed at higher redshifts, as is 
shown in Fig. [51 where the dimensionless ratios a /{He) 
and a /{He) with respect to the bare and dressed Hub- 
ble parameters are plotted. At last scattering, z ~ 1100, 
a/(iJc) -6 X 10^6. As t ^ 0, a/ {He) oc t^/^. 




FIG. 5: The magnitude of the dimensionless ratios a/{cH) 
(solid curve) and a/{cH) (dashed curve), where a{z), H{z), 
and H{z) are varied with redshift, for best-fit values of H^, 
.fvO- 

The object of estimating the relative deceleration scale 
was to determine whether its magnitude is physically ac- 
ceptable. This is the case. Although the physical magni- 
tude of the relative deceleration is larger at earlier epochs, 
over most of the history of the universe a is of the or- 
der of 10^^"ms^^, which is tiny. From Fig. [4l at z = 2 
we have a ~ 4.5 x 10~^''ms~^, which corresponds to an 
expansion age of 4Gyr, 27% of the current age of the uni- 
verse in wall time. By comparison, the Pioneer anomaly 
in the solar system [48| occurs at an acceleration scale of 
(8.74 ± 1.33) X 10^i°ms-2 a value which a attains 
only at a redshift of 3.85 equivalent to an expansion age of 
2.1Gyr in wall time. Furthermore, the relative decelera- 
tion scale should largely affect dynamics in the transition 
zones between walls and voids. At earlier epochs the void 
fraction is less: at z = 2 it is 44%, and at z = 3.85 it is 
28%. At z = 10, when a ~ 2 x lO'^ms-^, it is 10%. 

In the absence of an exact timelike symmetry of the 
background there is no obvious solution to the problem 
of how to keep two clocks synchronised in general relativ- 
ity. The CEP proposes a solution to this conundrum: the 
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evolving average density provides the relevant regional 
clock. Even though we are talking about weak fields in 
cosmology, and small relative decelerations between ex- 
panding regions of different densities, the fact that the 
relative decelerations are integrated over the lifetime of 
the universe means that the cumulative clock variance 
can be large for the density contrasts that are observed. 



B. Modified Newtonian dynamics? 

The fact that the present epoch value of a turns out 
to coincide with the MOND scale (46l . is intriguing. 
To start thinking about possible connections, it pays to 
recall that a is an estimate of the relative deceleration of 
wall regions, at finite infinity where space is still expand- 
ing, with respect to the cosmological volume average at 
any epoch. It is certainly an acceleration scale at which 
dynamical gradients in the Ricci curvature of space are 
likely to affect dynamics of particles between walls and 
void regions. 

Since we are no longer dealing with asymptotically 
flat geometries with exact timelike Killing vectors, it 
is quite possible that the solution to the Kepler prob- 
lem for bound geodesies in galaxies should be modified. 
Therefore the possibility that the MOND phenomenon 
is related is a reasonable hypothesis. However, the es- 
timate we have made of the relative deceleration scale, 
a, refers specifically to the finite infinity scale relative to 
the volume average. Given that the outskirts of galaxies 
are expected to lie within finite infinity, it means that 
no detailed quantitative comparisons can be made until 
the transition zone around finite infinity is more directly 
modelled. 

I will not attempt to look at the problem of the rotation 
curves of galaxies in the present paper, but will make a 
few observations that would follow if the effects of MOND 
are simply a modification of Newtonian dynamics in a 
static background which arise from dynamical density 
gradients affecting the Ricci curvature of space. 

The first important point is that, since we are dealing 
with an effect which is most pronounced between walls 
and voids, then the distance of a galaxy to a finite infinity 
boundary should play a role, at least naively. One might 
expect more pronounced effects for galaxies in filaments 
[5(j | as compared to rich clusters of galaxies in thick walls, 
where the distance to finite infinity is greater. In fact, 
rich clusters of galaxies often tend to harbour elliptical 
galaxies for which the MOND results differ little from 
Newtonian expectations. But apart from this, it appears 
that for the nearby redshifts over which it is tested, the 
MOND scale [iG^ is uniform for many different galaxy 
environments. This runs counter to naive intuition about 
the distance to finite infinity, and suggests that some key 
physical insight remains to be found. 

The second point is that, although the magnitude of 
the relative deceleration is larger at higher redshifts, 
the frequency of voids is also less. There may be a 



trade-off between these two competing factors as to 
an optimal redshift at which dynamical effects would 
be most evident. For the best-fit parameters [1^, the 
void fraction reaches 50% at a redshift z — 1.52 when 
a ~ 3.4 X 10~^°ms~^. It is possible that effects re- 
sulting from the relative deceleration of the background 
may be important for structure formation. Toy model 
calculations based on exact inhomogeneous solutions of 
Einstein's equations show that the nonlinear treatment 
of inhomogeneities in general relativity can considerably 
enhance the rate of structure formation 51 1 . 

Given that the expected distance between finite in- 
finity and the outskirts of a galaxy does not suggest a 
direct link between ao and amond, the close match be- 
tween ao and amond may be purely coincidental, even 
if both effects are related to Ricci curvature gradients 
which are typically of the same order of magnitude. As 
no direct link to MOND has yet been established, it is 
not possible to say whether the redshift dependence of a 
should be linked to a redshift dependence in the MOND 
scale. MOND is purely empirical and the fitting of ro- 
tation curves requires both an acceleration scale and an 
empirical function, fj.{x), which interpolates between the 
Newtonian and "modified dynamics" regimes. Any first 
principles treatment would need to describe the transi- 
tion zone between finite infinity and its interior, which 
might conceivably be related empirically to the interpo- 
lating function. We are at an early stage of gathering 
the pieces of a puzzle, where all observations have to be 
treated with rigorous scrutiny, and where no firm conclu- 
sions should be drawn until a compelling theoretical case 
is assembled. 

Finally, it should be noted that the best-fit parame- 
ters of Ref. indicate that nonbaryonic dark matter 
is likely to be the dominant component of matter in the 
universe by mass, even if the relative fraction is gener- 
ally somewhat lower than in the ACDM model. If there 
is any link to MOND, then it is the phenomenology of 
MOND that one might hope to explain, rather than an 
alternative to nonbaryonic dark matter. 



VI. WEYL CURVATURE AND INITIAL 
CONDITIONS 

The cosmological equivalence principle has been ap- 
plied here principally on the scale of macroscopic cosmo- 
logical averages. A question remains as to what extent 
we should take it as a universal principle? Given that the 
universe is well approximated by frames of the sort ([T]) at 
the earliest epochs which are observationally tested, such 
as the epochs relevant to big bang nucleosynthesis, it is 
certainly tempting to try to apply the CEP on arbitrar- 
ily small scales within the past light cone at the earliest 
epochs. For any perfect fiuid the pressure is determined 
by the density and no modification of the CEP is re- 
quired. However, if one is dealing with spinning fluids or 
gravitational waves which cannot be treated by a single 
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collective scalar degree of freedom, then modifications to 
the stated principle would be necessary. Is the process of 
"getting rid of local Ricci curvature" to recalibrate rods 
and clocks so that the effects of gravity disappear in a 
small region, as the equivalence principle demands, re- 
lated to the mathematical problem of Ricci flow? This 
is certainly a possibility, which has been discussed in the 
context of cosmological averaging by Buchert and Car- 
fora [HI. 

As far as the rest of the curvature degrees of freedom 
are concerned, the CEP in fact incorporates a strong 
statement about average Weyl curvature, since the GIF 
([1]) has vanishing Weyl curvature. Of course, Weyl curva- 
ture does not vanish everywhere identically - it is nonzero 
in the vicinity of collapsed structures such as stars and 
black holes, or in gravitational wave ripples. However, 
the statement of the CEP means that average effects of 
Weyl curvature need not be considered in the calibration 
of rods and clocks of generic cosmological frames. This 
"average Weyl curvature condition" may at first sight 
seem stronger than Penrose's "Weyl curvature hypothe- 
sis" jsSi] that the universe began in an initial state with 
total Weyl curvature exactly zero. However, at the epoch 
of last scattering the universe was very close to a stan- 
dard FLRW model, and all FLRW models have vanishing 
Weyl curvature. Thus the "average Weyl curvature con- 
dition" incorporated within the CEP could be viewed as 
a relatively weak phenomenological statement about the 
evolution of a universe whose initial state at last scatter- 
ing is consistent with the predictions of inflation. Alter- 
natively, it is also consistent with any other cosmological 
scenario which solves the flatness and horizon problems 
to give a close to spatially flat FLRW universe with near 
scale-free perturbations by the epoch of big bang nu- 
cleosynthesis, whether such a scenario satisfies Penrose's 
Weyl curvature hypothesis or not. 

It is my view, however, that in searching for potential 
scenarios for initial conditions, Penrose's Weyl curvature 
hypothesis needs to be taken seriously, and could be re- 
lated to a generalised cosmological equivalence principle, 
if it can be formulated to apply at earlier epochs in the 
very early universe. To understand my rationale let us 
recall that the Weyl curvature tensor includes any non- 
local curvature in a manifold, whereas the Ricci tensor 
encodes purely local curvature since it is directly related 
to the energy-momentum tensor via the Einstein equa- 
tions. 

Since general relativity is a local causal theory, given 
that the observable universe was in a global state very 
close to a FLRW geometry with zero Weyl curvature 
at last scattering, then the only Weyl curvature we are 
allowed today is that which has accumulated by local 
causal processes within the past light cone at any event: 
in particular, by gravitational collapse and production of 
gravitational waves. The Weyl tensor encodes tidal cur- 
vature information on local scales which grows as matter 
clumps. On the large scales where the universe is still ex- 
panding, Weyl curvature cannot be important in defining 



the average geometry, since the universe has only had a 
finite time over which to evolve from its state at last scat- 
tering. This is also the reason why there is a statistical 
scale of average homogeneity. 

If a version of the CEP can be taken to apply at even 
the earliest epochs, then it amounts to the statement 
that even throughout its earliest history the background 
universe contains no "nonlocal curvature" that cannot 
have evolved causally within the past light cone at any 
event. For as far back as a 4-dimensional spacetime con- 
tinuum has any meaning it would then make sense to 
be able to choose a CIF in which the average effects of 
density are volume-contracting. Large classes of models, 
such as Bianchi models with anisotropic flows, would be 
cosmologically irrelevant. If such a CEP should survive 
through the epoch of inflation or any other very early 
universe scenario, then it would coincide with Penrose's 
Weyl curvature hypothesis. Ultimately these conceptual 
issues might inform quantum gravity and quantum cos- 
mology. 



VII. DISCUSSION 

In this paper I have extended the strong equivalence 
principle to account for the average effect of the den- 
sity of matter in the dcflnition and relative calibration 
of clocks in inertial frames on cosmological scales. Since 
the resulting cosmological equivalence principle relates 
the single scalar degree of freedom of Newtonian gravity 
to the framework of general relativity, it may provide a 
means to better understand the calibration of cosmolog- 
ical weak fields once density perturbations have grown 
large to form a universe that is very inhomogeneous on 
scales of tens of megaparsecs. It should thereby give a 
setting for better understanding the Newtonian limit in 
the dynamical situation of cosmology. The numerical es- 
timate of the relative deceleration between observers in 
the walls around galaxy clusters and volume-average ob- 
servers in voids, typically of order lO^^^ms"^, is accept- 
ably small for weak field scales and yet leads cumulatively 
to the present epoch clock rate variance of 38% found in 
Refs. [Hill]. Intriguingly, at redshifts z < 0.25 the rela- 
tive deceleration required by the CEP coincides precisely 
with the empirical acceleration scale of MOND [46l . 

At a conceptual level I have attempted to present a 
framework for the consistent definition of average inertial 
frames in relation to average dynamically-varying mat- 
ter densities in cosmological general relativity. The hope 
is that the cosmological equivalence principle is thereby 
a key step to the incorporation of Mach's principle into 
general relativity, in the way that Einstein intended but 
never quite realised. Mach's principle is most commonly 
invoked in distinguishing inertial frames from rotating 
frames d, [H, [H, . What is studied here is a different 
aspect of Mach's principle: the role of the average volume 
deceleration of the local geometry in defining the stan- 
dard of time of inertial frames. In the absence of a time- 
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like Killing vector, the evolution of the average density 
provides a relevant clock. To fully incorporate Mach's 
principle in general relativity, it is of course necessary 
to deal with the other dynamical gravitational degrees of 
freedom which can affect the distinction of inertial frames 
from rotating ones, such as gravitational waves (56l |. 

In this paper I have expounded the view that to deal 
with the volume-contracting average dynamical effects 
of matter density, a reduction to a frame ([T]) is the rel- 
evant step in the normalisation of gravitational energy 
before the final step to a static Minkowski space. It is 
quite possible that when average degrees of freedom in 
addition to the scalar Ricci curvature are considered, in 
order to deal with gravitational waves and spinning mat- 
ter fluids, there are other steps in the relevant relative 
calibrations of rods and clocks. After all, energy, mo- 
menta, and angular momenta are only defined with re- 
spect to a frame. In the dynamical regime of general rel- 
ativity the question arises as to which collective average 
frames have physical utility in the absence of exact sym- 
metries described by Killing vectors. My view is that a 
truly deep understanding of quasilocal gravitational en- 
ergy and momentum is still to be found, but the path 
to such enlightenment requires a better conceptual un- 
derstanding of the equivalence principle in application to 
collective dynamical degrees of freedom of matter fields. 

Historically speaking, in the early stages of the devel- 
opment of general relativity Einstein did not fully appre- 
ciate the dynamical importance of the energy and mo- 
mentum of spacetinie itself. Spacetime is inevitably dy- 
namical for matter obeying the strong energy condition. 
Einstein's first journey through the conceptual landscape 
of cosmological general relativity had him worrying about 
boundary conditions at spatial infinity [lO| . as he over- 
looked the possibility that the universe had a beginning. 

Since general relativity is causal the geometry at any 
event can only depend on events within its past light 
cone, and is independent of what lies beyond the particle 
horizon. Thus boundary conditions at spatial infinity 



beyond the particle horizon are physically irrelevant if the 
universe had a beginning, a possibility that Einstein did 
not consider when he first formulated his static universe 
^ld\ . For a universe like ours which had a beginning, the 
initial conditions are of vital importance in determining 
the relevant weighted average of the apparent motions, 
as I have discussed in Sec. I VII 

The conceptual journey discussed in this paper arose 
in an effort to model the universe more realistically 
[ESU^i to account for the structure we actually observe, 
by realising that the quasilocal gravitational energy of a 
dynamical spacetime geometry - which has real effects on 
the calibration of clocks - should be an essential feature 
of a universe with large dynamical density gradients. If 
successful, this will eliminate the need for a cosmological 
constant or other fiuidlike vacuum energy as the source 
of "dark energy" , but it still leaves the other cosmologi- 
cal problem, why A = 0, unsolved. If we take the strong 
equivalence principle literally then A must be zero since 
otherwise we could not have a vacuum Minkowski space- 
time for our local inertial frames. My own personal view 
is that quantum field theoretic calculations based in a 
flat spacetime which suggest that A ~ Mpjg^j^^jj, miss the 
mark, because the spacetime vacuum cannot be under- 
stood without accounting for the intrinsically dynamic 
nature of spacetime. It is not a problem for flat space 
quantum fleld theory. While the cosmological constant 
problem is no doubt a problem for quantum gravity, I be- 
lieve that quantum gravity research might benefit from a 
more physical understanding of dynamical gravitational 
energy and the equivalence principle. 
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